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(54) Movement control in a lithographic apparatus 

(57) A method, system, and computer program 
product for non-real-time trajectory planning and real- 
time trajectory execution. A trajectory planning process 
receives data generated by high-level control software. 
This data defines positions and scan velocities, where 
multiple axis motion must be precisely synchronized. 
The trajectory planning process creates sequences of 
constant acceleration intervals that allow critical mo- 
tions to be executed at maximum throughput. The out- 
put of a trajectory planning process is known as a profile. 
A profile executor, using the profile output by the trajec- 
tory planner process, generates continuous, synchro- 
nized, filtered, multi-axis position and acceleration conrv 
mands (i.e., execution data) that drive control servos. 
Time intervals generated by the trajectory planner are 
quantized to be integer multiples of a real time clock pe- 
riod. The trajectory planner outputs have infinite jerk, but 
are smoothed by filters in the profile executor to both 
limit jerk and minimize servo-tracking errors. The trajec- 
tory planner allows time for the profile executor filters, 
but does not restrict fine tuning of the shape of these 
filters, provided that the width of the tuned filter does not 
exceed the allowed time. 



100 



Control process 



110 



J 



Trojectory planner 



Initial control 
data 120 130 



Trojectory data 
140 150 



Profile executor 



r 



Control system 



Execution dato 
160 170 



BEST AVAILABLE COPY 



FIG.1 



Printed by Jouve. 7S001 PARIS (FR) 



BNSOOCID: <EP. 



.1422569A2J_> 



EP 1 422 569 A2 



Description 

BACKGROUND OF THE INVErsTTION 

Field of the Invention s 

[0001] The invention described herein relates to tra- 
jectory planning and execution for movable components 
of a photolithography system. 

10 

Related Art 

[0002] Many Industrial processes involve motion 
along trajectories that are defined by precise positions 
at specific times. Photolithography is an example of i5 
such a process. In a photolithography process, an illu- 
mination source projects an Illumination beam. The 
beam passes through, or is reflected off of, a reticle to 
enable the transfer of a reticle image from the reticle to 
a substrate, such as a semiconductor wafer. 20 
[0003] Scanning techniques are employed In photoli- 
thography processes to project a reticle image into a 
substrate. These scanning techniques Involve moving 
a reticle across an illumination slot to allow the reticle 
image to be exposed onto a substrate that is simultane- 25 
ousty moving. Reticles and substrates are exposed on 
stages that are capable of motion in one or more dimen- 
sions. 

[0004] In general, digitally-calculated motion trajecto- 
ries enable high-precision multi-axis motion control sys- 30 
tems to achieve high throughput. Of special interest are 
the motion trajectories followed by the wafer and reticle 
in a step and scan lithography tool during the process 
in which the reticle pattern is imaged and exposed onto 
the wafer surface. The entire wafer surface is exposed 35 
in a sequence of field scans. The exposure of each field 
requires that the wafer and reticle be simultaneously 
scanned across the exposing field of the imaging optics 
at precisely synchronized, constant velocities. The ratio 
of the reticle and wafer velocities must exactly match 40 
the magnification of the imaging optics. 
[0005] To maximize system throughput, it is desirable 
to expose a wafer in the minimum possible time. In a 
typical trajectory used to expose a wafer, after each field 
exposure, the wafer stage must be asynchronously 45 
stepped from an initial state (i.e., a position and velocity) 
at the end of a field scan, to a new slate (i.e., a new 
position and velocity) at the start of the next field scan. 
Similariy, the reticle stage must also be asynchronously 
stepped from an initial state at the end of a field scan, so 
to a new state at the start of the next field scan. Conse- 
quently, minimizing the amount of time between scans 
would improve system throughput. A system and meth- 
od Is therefore needed to plan and execute optimal tra- 
jectories for photolithography components for the inter- ss 
vals between constant velocity scans. 



BRIEF SUMIWIARY OF THE INVENTION 

[0006] The invention described herein consists of a 
method, system, and computer program product for 
non-real-time trajectory planning and real-time trajecto- 
ry execution. A trajectory planning process receives da- 
ta generated by high-level control system. The high-lev- 
el control system can be implemented In software. This 
data defines positions and scan velocities, where mul- 
tiple axis motion must be precisely synchronized. The 
trajectory planning process creates sequences of con- 
stant acceleration intervals that allow critical motions to 
be executed at maximum throughput. The output of a 
trajectory planning process is known as a profile. A pro- 
file executor, using the profile output by the trajectory 
planner process, generates continuous, synchronized, 
filtered, multi-axis position and acceleration commands 
(I.e., execution data) that drive control servos. 
[0007] All time intervals generated by the trajectory 
planner are quantized to be integer multiples of a real 
time clock period. The trajectory planner outputs have 
infinite jerk, but are smoothed by filters in the profile ex- 
ecutor to both limit jerk and minimize servo-tracking er- 
rors. The trajectory planner allows time for the profile 
executor filters, but does not restrict fine tuning of the 
shape of these filters, provided that the width of the 
tuned filter does not exceed the allowed time. 
[0008] The combined trajectory planning and execu- 
tion processes allows the exposure sequence to be 
smoothly stopped and restarted after completing any in- 
dividual field scan. 

BRIEF DESCRIPTION OF THE FIGURES 

[0009] The foregoing and other features and advan- 
tages of the invention will be apparent from the follow- 
ing, more particular description of a preferred embodi- 
ment of the invention, as illustrated in the accompanying 
drawings. 

FIG. 1 is a block diagram illustrating the fundamen- 
tal modules of an embodiment of the invention, and 
the data flow between them. 
FIG. 2 illustrates the data flow between modules in 
greater detail, according to an embodiment of the 
invention. 

FIG. 3 illustrates a computing platfonn on whteh a 
software embodiment of the Invention can execute. 
FIG. 4 is a flowchart Illustrating, in general, the tra- 
jectory planning process, according to an embodi- 
ment of the invention. 

FIG. 5 is a flowchart illustrating the trajectory plan- 
ning process in greater detail, according to an em- 
bodiment of the Invention. 

FIG. 6 is a flowchart illustrating the process of cal- 
culating the motion of a component of a photolithog- 
raphy system for a given axis, according to an em- 
bodiment of the Invention. 
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FIG. 7 is a flowchart Illustrating the process of ad- 
justing a constant velocity scan length, according to 
an embodiment of the Invention. 
FIG. 8 is a flowchart illustrating the process of cal- 
culating an acceleration interval, taking into account s 
maximum velocity and acceleration constraints, ac- 
cording to an embodiment of the invention. 
FIG. 9 illustrates a profile that contains two accel- 
eration intervals and has no maximum velocity con- 
straint. 10 
FIG. 1 0 Is a flowchart Illustrating the process of cal- 
culating acceleration intervals in the case where 
there is no maximum velocity constraint, according 
to an embodiment of the invention. 
FIG. 11 Illustrates a profile that contains three ac- is 
celeration Intervals and having a maximum velocity 
constraint. 

FIG. 12 Is a flowchart Illustrating the process of cal- 
culating acceleration Intervals in the case where 
there is a maximum velocity constraint, according ^ 
to an embodiment of the invention. 
FIG. 13 is a flowchart illustrating the process of cre- 
ating a zero velocity point between scans, accord- 
ing to an embodiment of the Invention. 
FIG. 14 is a flowchart Illustrating the process of ad- 
justing zero velocity time Intervals to synchronize 
multiple scan axes, according to an embodiment of 
the Invention. 

FIG. 1 5 is a flowchart Illustrating the process of cre- 
ating a delay at the start of a motion trajectory, to 30 
compensate for fitter delays, according to an em- 
bodiment of the invention. 

FIG. 1 6 is a flowchart illustrating the process of pro- 
file execution, according to an embodiment of the 
invention. 35 
FIG. 1 7 is a flowchart illustrating the process of fil- 
tering an integrated position signal, according to an 
embodiment of the invention. 
FIG. 18 is a flowchart Illustrating the process of in- 
terpolating a filtered position signal, according to an ^ 
embodiment of the invention. 
FIG. 19 is a flowchart illustrating the process of fil- 
tering an up-sampled signal to create an interpolat- 
ed signal, according to an embodiment of the inven- 
tion. 45 
FIG. 20 is a flowchart illustrating the process of In- 
terpolating an acceleration state signal, according 
to an embodiment of the Invention. 
FIG. 21 is a block diagram illustrating a profile ex- 
ecutor and its components, according to an embod- so 
iment of the invention. 

FIG. 22 is a block diagram Illustrating a filter for an 
integrated position signal, according to an embodi- 
ment of the Invention . 

FIG. 23 is a block diagram illustrating an interpola- S5 
tion component of a profile executor, according to 
an embodiment of the invention. 



DETAILED DESCRIPTION OF THE INVENTION 

[0010] A prefen-ed embodiment of the present inven- 
tion is now described with reference to the figures where 
like reference numbers indicate identical or functionally 
similar elements. Also in the figures, the left most digit 
of each reference number corresponds to the figure in 
which the reference number is first used. While specific 
configurations and arrangements are discussed, it 
should be understood that this is done for Illustrative pur- 
poses only. A person skilled In the relevant art will rec- 
ognize that other configurations and arrangements can 
be used without departing from the spirit and scope of 
the Invention. It will be apparent to a person skilled in 
the relevant art that this invention can also be employed 
In a variety of other devices and applications. 

I. Overview 

[0011] The invention described herein comprises a 
non-real-time planning algorithm, and a real time exe- 
cution algorithm. The trajectory planning process re- 
ceives data generated by high-level control software. 
This data defines positions and scan velocities, where 
multiple axis motion must be precisely synchronized. 
The trajectory planning process creates sequences of 
constant acceleration Intervals that allow critical scan- 
ning motions to be executed at maximum throughput. 
The profile executor, using trajectory planner outputs (i. 
e., the profile), generates continuous, synchronized, fil- 
tered, multi-axis position and acceleration commands 
that drive control servos. All time intervals generated by 
the trajectory planner are quantized to be integer multi- 
ples of a real time clock period. The trajectory planner 
outputs have infinite jerk, which Is smoothed by the fil- 
ters in the profile executorto both limit jerk and minimize 
servo-tracking errors. The trajectory planner allows time 
for the profile executor filters, but does not restrict fine 
tuning of the shape of these filters, provided that the 
width of the tuned filter does not exceed the allowed 
time. 

[0012] The invention is illustrated generally in FIG. 1 . 
A control process 110 Is responsible for providing initial 
control data 120 to a trajectory planner 130. Initial con- 
trol data 120 includes Information necessary to begin 
trajectory planning for an interval between constant ve- 
locity scans. Because this interval connects two scans, 
it is sometimes referred to as a link. Initial control data 
120 includes the starting point of a scan, the required 
distances and velocities of a scan, and stopping points. 
Given the initial control data 120, trajectory planner 130 
then produces trajectory data 140, also known as a pro- 
file. Trajectory data 140 includes definition of quantized 
time Intervals of constant acceleration for the link, as 
well as the acceleration during these intervals. Trajec- 
tory data 140 also includes the position and velocity at 
the start of each acceleration interval. Trajectory data 
140 is then used by profile executor 150. Profile execu- 
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tor 150 uses trajectory data 140 to produce execution 
data 160. Execution data 160 is a real time data streann 
defining positions and accelerations. Execution data 
1 60 represents a series of commands to control system 
1 70. Control system 1 70 represents the electromechan- 
ical control mechanism used to move components of the 
photolithography system (e.g., the wafer stage and ret- 
icle stage) according to execution data 160. 
[0013] FIG. 2 is a more detailed block diagram of the 
system of the invention. System 200 includes trajectory 
planner 130, profile executor 150, and control system 
170. Trajectory planner 130 designs device trajectories 
in accordance with control data 120. This design results 
in the generation of trajectory data 140. Trajectory data 
140 can be viewed as a state signal set that defines a 
trajectory for a device, such as a trajectory for a sub- 
strate stage or a reticle stage in a photolithography tool. 
Planner 130 can be Implemented as a computer pro- 
gram product that is executed by a computer. However, 
planner 130 can also be implemented in hardware and/ 
or firmware or a combination thereof. In an embodiment 
of the invention, trajectory planner 130 designs these 
trajectories in an off-line process. Thus, In the case of 
photolithography, trajectory data 140 is generated prior 
to commencement of substrate stage or reticle stage 
scanning operations. 

[0014] Trajectory data 140 defines a plurality of con- 
stant acceleration states that govern the motion of a de- 
vice, such as a substrate stage or a reticle stage, that 
is included in control system 170. Trajectory data 140 
Includes a position state signal 21 0, a velocity state sig- 
nal 212, and an acceleration state signal 214. For each 
constant acceleration state, state signals 21 0, 212, and 
214 have corresponding values. These constant accel- 
eration states exist for variable length time intervals that 
are integer multiples of a slow clock period such as 1 
millisecond. 

[0015] As described above, trajectory data 140 in- 
cludes position state signal 110, velocity state signal 
112, and acceleration state signal 114. These state sig- 
nals generated by trajectory planner 140 define motion 
in a single dimension for a device, such as a substrate 
stage or a reticle stage. 

[0016] Relationships between these signals can be 
expressed as P = Pf-^"^ + and V= + At. In 
these expressions, where Pf represents position state 
signal 21 0, P represents a position signal, A represents 
acceleration state signal 214, Vf represents velocity 
state signal 212, and represents a velocity signal, Po- 
sition state signal 210 is a non-equidistantly sampled 
second order curve. Accordingly, velocity state signal 
212 is a non-equidistantly sampled trapezoidal signal. 
[001 7] Acceleration state signal 21 4 includes a series 
of square and/or rectangular pulses that can have any 
amplitude and pulse width. These pulses correspond to 
transitions in velocity state signal 212 values. At the ris- 
ing and falling edges of these pulses, acceleration state 
signal 214 has a first derivative that is large. This first 



derivative is referred to herein as "jerk". Since these ris- 
ing and falling edges provide for square and/or rectan- 
gular pulses, acceleration state signal 214 has jerk val- 
ues that are infinity (i.e., delta functions). 
s [0018] The values of velocity state signal 212 and ac- 
celeration state signal 21 4 only change at specific times 
deflnedbytrajectory planner 130. Thus, instead of being 
represented by sets of samples that are unifomily 
spaced In time, these signals are represented by values 
10 corresponding to specific moments in time when accel- 
eration state signal 214 changes. As described above, 
these moments are separated by time intervals that are 
integer multiples of the slow clock period. 
[0019] Referring again to FIG, 1 , profile executor 1 50 
IS receives trajectory data 1 40 from trajectory planner 1 30 
and converts these signals Into an output signal set, ex- 
ecution data 160. Execution data 160 includes a real 
time position signal and a real time acceleration output 
signal. These output signals correspond to trajectory da- 
20 ta 1 40, but are discrete time signals occurring at a fixed 
fast clock rate. An exemplary fast clock rate is 1 6 kHz. 
Other clock rates can also be used. Discrete time sig- 
nals occurring at this fast clock rate have values spaced 
at 0.0625 millisecond time intervals. The process of in- 
25 creasing the clock rate of a discrete time signal is re- 
ferred to herein as interpolation. Profile executor 150 
sends execution data 1 60 to control system 1 70. 
[0020] In generating execution data 1 60, profile exec- 
utor 1 50 filters the state signals of trajectory data 1 40 to 
30 reduce excessive jeric, such as those described above. 
In so doing, profile executor 160, through filtering, 
avoids exciting certain resonances In control system 
1 70 and produces discrete time signals at the fast clock 
rate. Since this fast clock rate Is greater than the slow 
35. clock rate associated with trajectory data 1 40, profile ex- 
ecutor 150 perfomns filtering and Interpolation of state 
signals to produce execution data 160. 
[0021] Execution data 160 includes a position output 
signal 21 8 and an acceleration output signals 220. Like 
40 state signals 210-214, these signals are discrete sig- 
nals. However, unlike state signals 210-214, position 
signals 21 8 and acceleration signal 220 each have val- 
ues that are spaced at uniform time intervals. These unl- 
fomri time intervals are equal to one period of the fast 
45 clock rate described above. 

[0022] Profile executor 1 50 Is Implemented as a com- 
puter program product that can be executed by a com- 
puter system. In an embodiment of the invention, this 
computer is a fixed precision arithmetic micro-controller, 
50 such as a SHARC® ADSP-21 062 produced by Analog 
Devices, Inc. of Norwood, Massachusetts. 
[0023] The functions perfonned by trajectory planner 
130, whk^h is capable of calculating trajectories as 
shown in process 500, are preferably implemented in 
55 software. Alternatively, the same may be implemented 
using hardware or a combination of hardware and soft- 
ware. 

[0024] An example of a computer system 300 is 
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shown generally in FIG. 3. Such a computer system can 
be used to perform the processing of trajectory planner 
130 or profile executor 150. The computer system 300 
represents any single-processor computer. Single- 
threaded and multi-threaded computers can be used, s 
Unified or distributed memory systems can be used. 
[0025] The computer system 300 includes one or 
more processors, such as processor 304. One or more 
processors 304 can execute software Implementing the 
operations described below in FIG. 5. Each processor io 
304 is connected to a communication bus 302 (e.g., 
cross-bar or network). Various software embodiments 
are described in terms of this exemplary computer sys- 
tem. After reading this description, it will become appar- 
ent to a person skilled in the relevant art how to imple- is 
ment the Invention using other computer systems and/ 
or computer architectures. 

[0026] Computer system 300 also includes a main 
memory 306, preferably random access memory 
(RAM), and can also include a secondary memory 308. so . 
The secondary memory 308 can include, for example, 
a hard disk drive 310 and/or a removable storage drive 
312, representing a floppy disk drive, a magnetic tape 
drive, an optical disk drive, etc. The removable storage 
drive 312 reads from and/or writes to a removable stor- 25 
age unit 314 in a well known manner. Removable stor- 
age unit 314 represents a floppy disk, magnetic tape, 
optical disk, etc., which is read by and written to by re- 
movable storage drive 312. As will be appreciated, the 
removable storage unit 31 4 includes a computer usable 30 
storage medium having stored therein computer soft- 
ware and/or data. 

[0027] In alternative embodiments, secondary mem- 
ory 308 can include other means for allowing computer 
programs or other instructions to be loaded into compu- 35 
ter system 300. Such means can include, for example, 
a removable storage unit 322 and an interface 320. Ex- 
amples can include a program cartridge and cartridge 
interface (such as that found in video game devbes), a 
removable memory chip (such as an EPROM, or 40 
PROM) and associated socket, and other removable 
storage units 322 and interfaces 320 which allow soft- 
ware and data to be transferred from the removable stor- 
age unit 322 to computer system 300. 
[0028] Computer system 300 can also Include a com- 45 
municatlons interface 324. Communications Interface 
324 allows software and data to be transferred between 
computer system 300 and external devices via commu- 
nications path 326. Examples of commun cations inter- 
face 320 can include a modem, a network interface so 
(such as Ethemet card), a communications port, etc. 
Software and data transferred via communications in- 
terface 324 are in the form of signals which can be elec- 
tronic, electromagnetic, optical or other signals capable 
of being received by communications interface 324, via 55 
communications path 326. Note that communications 
interface 324 provides a means by which computer sys- 
tem 300 can interface to a network. 



[0029] The present invention can be implemented us- 
ing software running (that is, executing) in an environ- 
ment similar to that described above with respect to FIG. 
3. In this document, the temn "computer program prod- 
uct" is used to generally refer to removable storage unit 
314, a hard disk installed in hard disk drive 310, or a 
carrier wave carrying software over a communteation 
path 326 (wireless link or cable) to communication in- 
terface 324. A computer useable medium can include 
magnetic media, optical media, or other recordable me- 
dia, or media that transmits a carrier wave. These com- 
puter program products are means for providing soft- 
ware to computer system 300. 

[0030] Computer programs (also called computer 
control logic) are stored in main memory 306 and/or sec- 
ondary memory 308. Computer programs can also be 
received via communications Interface 324. Such com- 
puter programs, when executed, enable the computer 
system 300 to perform the features of the present inven- 
tion as discussed herein. In particular, the computer pro- 
grams, when executed, , enable the processor 304 to 
perform the features of the present invention. Accord- 
ingly, such computer programs represent controllers of 
the computer system 300. 

[0031] In an embodiment where the invention is im- 
plemented using software, the software can be stored 
In a computer program product and loaded into compu- 
ter system 300 using removable storage drive 312, hard 
drive 31 0, or communications Interface 324. Alternative- 
ly, the computer program product can be downloaded 
to computer system 300 over communications path 324. 
The control logic (software), when executed by the one 
or more processors 304, causes the processor(s) 304 
to perform the functions of the invention as described 
herein. 

[0032] In another embodiment, trajectory planner 130 
and/or profile executor 150 are implemented primarily 
in fimiware and/or hardware using, for example, hard- 
ware components such as application specific integrat- 
ed circuits (ASICs). Implementation of a hardware state 
machine so as to perform the functions described herein 
will be apparent to persons skilled in the relevant art(s). 

II. Trajectory Planning 

[0033] FIG. 4 shows a flowchart depicting an embod- 
iment of the operation of the trajectory planning process 
of the present invention. FIG. 4 shows how a trajectory 
is calculated, generally, by trajectory planner 130. The 
process begins with step 402. In step 404, a trajectory 
planner calculates a trajectory subject to velocity and 
acceleration constraints. These constraints represent 
physical limitations of control system 170. The calcula- 
tion of step 404 is described in greater detail below. In 
step 406, the trajectory calculated in step 404 is quan- 
tized. This represents decomposition of the trajectory in- 
to discrete time intervals. The process concludes at step 
408. 
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[0034] Trajectory planning is illustrated in greater de- 
tail in FIG. 5. The process begins with step 510. In step 
620, the trajectory planner receives Initial control data. 
In step 530, the preliminary trajectory of a given compo- 
nent of the photolithography system is calculated for 
each axis. At step 540, adjustments are made to the tra- 
jectory at the time Intervals in which velocity is zero, so 
as to synchronize motion in the different scan axes. In 
step 550, delays are inserted at the start of motion tra- 
jectories to compensate for filter delays that may occur 
in the profile executor 150. The process concludes at 
step 560. 

A. Calculation of Motion in Each Axis 

[0035] The process of calculating the motion of a com- 
ponent of a photolighography system In each axis, step 
530, Is illustrated In FIG. 6. The process begins with a 
step 61 0. In step 620, the length of any constant velocity 
scan is adjusted to account for the width of filters in the 
profile executor, settling time, and quantization. Specif- 
ically, adjustment of a constant velocity scan length is 
performed by extending end points of the constant ve- 
locity scan. "Settling time" refers to time required for ser- 
vo errors to reach an acceptable performance level. 
When a motion profile consists of acceleration intervals 
during which servo errors are not critical, followed by a 
time period of zero or constant velocity where servo er- 
rors must be minimized, it Is often necessary to allow a 
time period, after the acceleration is over, for settling. 
This settling time allows the servo errors to reach the 
appropriate performance level. Step 620 will be de- 
scribed in greater detail below. 
[0036] In step 630, acceleration intervals are calculat- 
ed for moving a component from position and velocity 
states at the end of one scan to the position and velocity 
state at the beginning of the next scan. It is desirable to 
do this in a minimum amou nt of time consistent with con- 
trol system constraints for maximum velocity and accel- 
eration. 

[0037] In step 640, a zero velocity point is identified, 
or if necessary created in the motion profile to provide 
a point where time can be added to the profile for the 
purpose of either synchronizing multiple axes, or to al- 
low the trajectory to be easily stopped and restarted af- 
ter any individual field scan. 

[0038] In step 650, acceleration intervals are quan- 
tized to produce a series of discrete acceleration inter- 
vals. Here, constant velocity time Intervals are rounded 
to an integer number of clock periods, and new accel- 
eration values are calculated to produce the same net 
motion that would have resulted If the non-integer ac- 
celeration periods had been used. This quantization 
procedure must comply with constraints for maximum 
velocity and acceleration, and, as a result, the total 
quantized time intervals may be longer, but never short- 
er than the Ideal, non-quantized intervals. The motion 
calculation process concludes at step 660. 



[0039] Step 620 above, the adjustment of constant 
velocity scan length, is illustrated in greater detail in FIG. 
7. The process begins at step 71 0. At this point, the start 
and end points of a constant velocity scan, as provided 

s in the initial control data, have been received. In step 
730, a determination Is made as to whether the respec- 
tive widths of all filters in the profile executor are accom- 
modated. If not, in step 740 the length of the constant 
velocity scan is extended to accommodate the widths of 

10 the filters. Not fully accommodating the filters is an op- 
tion that will result In part of the exposure scan not being 
at constant velocity. If both reticle and wafer filters are 
matched then synchronization of the stages during ex- 
posure can still be maintained. A varying exposure ve- 

'5 locity will require a dynamic correction of exposing pow- 
er to maintain dose unifonmity. In step 750, a determi- 
nation is made as to whether servo settling is accom- 
modated. If not, in step 760, the end points of the scan 
lengths are extended to accommodate servo settling. In 

20 step 770, a determination is made as to whether the start 
and end points are separated by a fixed number of pe- 
riods. If not, the process continues at step 780. At this 
step, start and end points are moved, i.e. scan lengths 
are extended, until separation between the points is a 

25 fixed number of periods. In step 790, the process con- 
cludes. 

[0040] Step 630, the computatiori of acceleration in- 
tervals, is illustrated in greater detail in FIG. 8. The proc- 
ess begins with step 810. In step 820, the end points of 

30 the scans, as adjusted in step 620, are received. In step 
830, optimum intervals of maximum acceleration are 
calculated to allow a component to move from the posi- 
tion and velocity state at the end of one scan to the po- 
sition and velocity state at the beginning of the next 

35 scan. A maximum velocity constraint is not yet applied. 
Step 830 is illustrated in greater detail in FIGS. 9 and 
10 below. 

[0041 ] In step 840, the two-state profile is modified in- 
to a three-state profile that is consistent with a maximum 

40 velocity (Vn^g^) constraint. FIGs. 1 1 and 1 2 illustrate this 
process, and will be detailed below. In the event the ve- 
locity exceeds the maximum value feasible for the con- 
trol system, a zero acceleration Interval is added be- 
tween the two maximum acceleration intervals, and the 

^5 lengths of the maximum acceleration intervals are mod- 
ified to maintain the same start and end points. The 
process concludes at step 860. 

[0042] An embodiment of step 830 is illustrated in 
greater detail in FIGS. 9 and 1 0. FIG. 9 shows a "bang- 
so bang" profile with two constant acceleration intervals 
that are consistent with a maximum acceleration <a,|,^) 
constraint. Graph 900 shows acceleration versus time. 
Initially, acceleration is zero, while position and velocity 
are given by and respectively. Motion is then ac- 
55 celerated to a level denoted with the variable a. Accel- 
eration a continues for a duration t^, before changing to 
acceleration -a. Acceleration -a continues for a duration 
t2. Finally, acceleration returns to zero, at which point 
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position and velocity are given by P2 and V2 respective- 
ly. 

[0043] FIG. 10 is a flowchart illustrating the process 
of calculating an acceleration tinne interval without a 
maximum velocity constraint, according to an embodi- 
ment of the invention. The process starts with step 1 006, 
in which V^.. Pi, V2, P2, and a^a^ are input. In step 1 01 0, 
the total displacement d is calculated, d = P2 - Pi- In 
step 1015, an acceleration variable a Is set to a^ax* 
[0044] In step 1 020, the quadratic equation is applied 
to detennine one of the two possible solutions for t^, 
shown here as t^^. The calculation Is as follows: 



'1.= 



-V,-t-[(VaH V;')/2+aj]'" 



10 



15 



*2e . 'ir ^ 



[0047] 



In step 1 040, the second solution is calculated, 



-V,-[(v/ + V;^)/24a^f 



The corresponding value for t2, shown here as t2d> is 
calculated on the basis of t^d* 



The corresponding value for t2, shown In step 1020 as 
t2ai is calculated on the basis of t^^: 



20 



[0045] In step 1025, the second solution for t., is cal- 
culated, identified as t.||, : 



The corresponding value for t2, shown here as t2b. is 
calculated on the basis of X^^: 



[0046] In step 1030, the variable a is set to -a^ax- A® 
In step 1020, the quadratic equation Is applied in step 
1 035 to determine one of the two possible solutions for 
t.|, shown here as t^i^. The calculation is as follows: 



-V, + [(Va' + V,')/2+arf]'" 



25 



30 



35 



40 



45 



SO 



[0048] In step 1 045, the number of acceleration states 
(N) is set to 2. In step 1050, one of the above t.| values 
is chosen, along with its corresponding t2 value. The se- 
lection is performed such that both t.| and t2 are positive 
real numbers, and the sum t^ + 12 is minimized. In step 
1055, t.|, t2, a, and N are returned. 
[0049] An embodiment of step 840, the application of 
velocity constraint to the computation of an accel- 
eration interval, is illustrated in greater detail in FIGS. 
11 and 12. Graph 1100 of FIG. 11 shows how accelera- 
tion changes over time, given a maximum velocity con- 
straint. Initially, acceleration is zero, and position and ve- 
locity are given by P^ and respectively. Acceleration 
then begins, quantified as acceleration a, and lasts for 
a duration Acceleration then returns to zero for a 
duration t^j^i.' after which position and velocity are given 
by P^y and V^y respectively. Another acceleration In- 
terval then occurs, where the acceleration has a value 
-a, for duration t2a- After t2a, the position and velocity are 
given by P2 and V2 respectively. 
[0050] FIG. 12 illustrates the calculation of accelera- 
tion intervals given a constraint V^ax> according to an 
embodiment of the invention. The process begins with 
step 1205, in which values for P^, V.,, P2, V2, a, t^, t2, 
and V^ax are input. In step 1210, V^^ is calculated to 
be V.| + at^. In step 1215, a determination Is made as to 
whether the absolute value of V^y exceeds V^qj^. If not, 
then the velocity constraint is not violated, and the proc- 
ess concludes at step 1220. If V,,,^ exceeds V^q^* then 
in step 1 225 a determination is made as to whether V^i^ 
exceeds zero. If so, then \fj^^ is setto M^^^ in step 1230. 
Otherwise, V„,^ Is set to -M^^^ in step 1235. 
[0051] In step 1 240, t.13 Is calculated, 



55 



The corresponding value for t2, shown here as t2g, is 
calculated on the basis of t^^: 
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In step 1245, ^ is calculated, 



[0052] In step 1250, the displacement d between po- 
sitions Pi and P2 Is calculated. In step 1255. the dis- 
tance traveled during t^^ Is calculated, 



2 



In step 1260, the distance traveled during tg^ is calcu- 
lated, 



2 



10 



15 



20 



In step 1265, the distance traveled during t^^ is calcu 
tated, 



[0053] In step 1270, P^,^ Is calculated, 



I- 25 



In step 1275, the duration X^i^^^ is calculated, 



30 



35 



is brolcen into two less optimal segments, wherein the 
first ends with velocity equal to zero, and the second 
starts with velocity equal to zero. The process concludes 
at step 1360. 

B. Adjusting Zero Velocity Time Intervals to 
Synchronize Multiple Scan Axes 

[0055] Step 540, the step of adjusting zero velocity 
time Intervals to synchronize multiple axes, is described 
in greater detail in FIG. 1 4. At this point in the trajectory 
planning process, all the required motion trajectories 
have been reduced to quantized intervals of constant 
acceleration, and furthermore there are defined zero ve- 
locity points between all synchronized portions of the 
trajectories. Using Integer arithmetic, additional time 
can be Inserted during the non-synchronized portions of 
the motion trajectories, at these zero velocity points. 
This allows the axes with shorter trajectory durations (i. 
e., non-critical paths) to be delayed to take the same 
time as the longest (i.e., critical) path axis. This causes 
the scanning motions to be precisely synchronized, with 
no loss in throughput. The process for each non-critical 
axis trajectory begins with step 1410. in step 1420, a 
non-critical axis trajectory is identified. In step 1430, the 
difference in time between the critical axis trajectory and 
the non-critical axis trajectory is calculated. In step 
1440, a delay is introduced, wherein the delay is equal 
to the difference, calculated in step 1430, at the zero 
velocity point of the non-critical axis trajectory. The proc- 
ess concludes at step 1450. 

[0056] Note that in alternative embodiments of the in- 
vention, the magnitude of the acceleration of an axis 
having a shorter trajectory duration can be reduced. 
Such reduction in acceleration can also be used in syn- 
chronizing multiple scan axes. 



40 



In step 1280, the number of acceleration intervals N is 
set to three. The process concludes at step 1 285 with 
the return of t^a, t^w. ^23. Pmid. and N, 
[0054] Step 640, the step of defining zero velocity 45 
points between scans, is illustrated in greater detail in 
FIG. 1 3. The process begins at step 131 0. in step 1 320, 
a determination is made as to whether the trajectory link- 
ing the end of one scan with the start of another scan 
starts or ends with zero velocity. If so, the process con- so 
eludes at step 1 360. if not, the process continues to step 
1 325, where a determination is made as to whether this 
trajectory has a reversal point. If so, then processing 
continues at step 1330. Here, the time of reversal is cal- 
culated. In step 1 340, this trajectory is split into two seg- 55 
ments at the reversal point (at which velocity is zero). 
Altematively, if the trajectory has no reversal point, then 
the process continues at step 1350. Here, the trajectory 



C. Compensation For Executor Fitter Delays 

[0057] In an embodiment of the invention, an addition- 
al step in the planning process is to account for any de- 
lays that may be introduced by the filtering used in the 
profile executor. In general, any time a new motion tra- 
jectory is started, axes with shorter filter delays can be 
intentionally delayed by an integer number of clock pe- 
riods. This allows synchronization of the mutti-axis mo- 
tion profiles produced by the profile executor. This proc- 
ess of compensating for filter delays is illustrated In 
greater detail in FIG. 15 and corresponds to step 550 of 
FIG. 5, The process begins with step 1510. Instep 1520, 
the axis trajectory having the longest execution filter de- 
lay Is identified, in step 1530, an axis trajectory having 
a shorter execution filter delay is identified. Instep 1540, 
the difference between the delays is calculated, in step 
1 550, a delay, equal to this difference, is introduced into 
the axis trajectory with the shorter execution filter delay. 
The process concludes at step 1560. 
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II. Profile Execution 

[0058] The profile executor integrates the output of 
the trajectory planner Into an infinite jerk position profile, 
and then uses filters to optimize the shapes of the posi- s 
tion profiles. In an embodiment of the invention, finite 
impulse response (Fl R) filters are used. Alternatively, in- 
finite impulse response (IIR) filters are used. Separate, 
independentiy-tuned filters can be used to optimize the 
shape of acceleration commands. The filters used to io 
smooth the infinite jerk position profiles must have 
matched group delays on all axes, unless the trajectory 
planner process accounts for the differences. Much of 
the executor process can be executed at a relatively 
slow clock rate, and then up-sampled to a higher clock is 
rate using an interpolating algorithm. The outputs of the 
profile executor are then fed to the control system. An 
Implementation of a system and method for profile exe- 
cution is described below and in U.S. Patent Application 
No. 09/757,622, "Method and System for Efficient and 
Accurate Filtering and Interpolation." filed January 11, 
2001 , and Incorporated herein by reference in its entire- 
ty. 

[0059] An embodiment of the profile execution proc- 
ess is illustrated in FIG. 16. The process begins with 
step 1 610. In step 1 620, trajectory data is received from 
the trajectory planner. As described above, the trajecto- 
ry data can be viewed as a set of state signals. In step 
1 630, an integrated position signal is created, based on 
the trajectory data. The position P^ at the time interval so 
n is given by the equation 

P^ = P| + nT (v, + nTA/2), 

35 

where 

P, is the initial position provided by signal 210, 
n is the sample time, 
T is the clock period, 

Vj Is the initial velocity given by signal 212, and ^0 
A is the acceleration, as given by signal 214. 
In step 1640, the integrated position signal is filtered, 
given an initial position signal, to create a filtered posi- 
tion signal. This step Is described in greater detail below. 
In step 1650, the filtered position signal Is interpolated, 
to create an interpolated signal, given an initial position 
signal. This step is likewise described in greater detail 
below. In step 1660, the interpolated signal is delayed 
to compensate for propagation delay and phase shift. 
The process concludes at step 1 670. so 
[0060] Step 1640, the filtering of an integrated posi- 
tion signal, is described in greater detail in FIG. 1 7. The 
process begins at step 1710. In step 1720, an absolute 
position array is applied to the integrated position signal, 
to create a parallel position signal set. in step 1 730, the 55 
initial position, P|. is subtracted from each parallel posi- 
tion signal of the created parallel position signal set. This 
creates a delta parallel position signal set. In step 1 740, 



the delta parallel position signal set is multiplied by filter 
coefficients, to create a parallel product signal set. In 
step 1 750. each of the signals of the parallel product 
signal set is summed to create a single scalar delta fil- 
tered signal. In step 1 760, a delayed initial position sig- 
nal T| is added to the delta filtered signal, to create the 
filtered position signal. The process concludes at step 
1770. 

[0061] Step 1650, the interpolation of a filtered posi- 
tion signal, is illustrated in greater detail in FIG. 1 8. The 
process begins at step 1810. In step 1820, the filtered 
position signal is upsampted. In step 1830, the upsam- 
pled signal is filtered to create the interpolated signal. 
[0062] In parallel with steps 1630 through 1660, ac- 
celeration state signal 214 Is interpolated separately. 
This takes place in step 1 663. In this step, acceleration 
state signal 214 is first up sampled at an interpolation 
rate L. The unsampled acceleration signal is then fil- 
tered. Filtering prevents high frequency resonance com- 
ponents from being introduced by the upsampling proc- 
ess. The resulting interpolated acceleration signal Is 
then delayed in step 1 666. The delay represents adjust- 
ment of the interpolated acceleration signal to address 
propagation delay and phase alignment. The process 
concludes at step 1670. 

[0063] Step 1830 above, the filtering of an up-sam- 
pled signal, is illustrated in greater detail in FIG. 1 9. The 
process begins at step 1910. In step 1920, an absolute 
position array is applied to the up-sampled signal to cre- 
ate a parallel up-sampled position signal set. In step 
1930, the initial position signal, Pj, is subtracted from 
each parallel up-sampled position signal, to create a del- 
ta parallel up-sampled position signal set. In step 1940, 
the delta parallel position signal set is multiplied by filter 
coefficients, to create a parallel product up-sampled sig- 
nal set. In step 1950, each of the signals of the parallel 
product signal set is summed, to create a single scalar 
delta filtered up-sampled signal. In step 1960, the de- 
layed initial position signal, P|, is added to the delta fil- 
tered up-sampled signal, to create an interpolated sig- 
nal. The process concludes at step 1970. 
[0064] Step 1 663 of FIG. 16, the step of interpolating 
the acceleration state signal, is described in greater de- 
tail in FIG. 20. The process begins with step 2010. In 
step 2020, the acceleration state signal Is upsampled. 
in step 2030, the upsampted signal is filtered to create 
an interpolated acceleratton signal. The process con- 
cludes at step 2040. 

[0065] Profile executor 1 50 can be Implemented as a 
computer program product that is executed by a com- 
puter. For example, the logic described above can be 
executed on a fixed precision arithmetic micro-control- 
ler, such as a SHARCXS) ADSP-21062 produced by An- 
alog Devices, Inc. of Norwood, Massachusetts. 
[0066] FIG. 21 is a block diagram of an implementa- 
tion of profile executor 150. Executor 150 includes a 
state signal interface 2102, an output signal Interface 
21 04, a slow clock portion 21 06, and a fast clock portion 
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21 08. State signal interface 2 1 02 receives trajectory da- 
ta 140 from trajectory planner 1 30. As described above 
with respect to FIG. 2, trajectory data 1 40 comprises po- 
sition state signal 21 0, velocity state signal 212, and ac- 
celeration state signal 21 4 Output signal interface 21 04 5 
sends execution data 160 to control system 170. 
[0067] Spanning slow and fast docl< portions 2106 
and 2108 are an acceleration signal patli 2112 and a 
position signal path 21 1 4. Acceleration signal path 21 1 2 
receives acceleration state signal 21 4 from state signal io 
interface 21 02 and generates acceleration output signal 
220. Position signal path 2114 receives state signals 
210, 212, and 214 (i.e., trajectory data 140) from state 
signal interface 21 02. Based on these received signals, 
position signal path 2114 generates position output sig- t5 
nal218. . 
[0068] Both acceleration signal path 2112 and posi- 
tion signal path 2114 include a plurality of signal 
processing elements. These elements perform opera- 
tions, such as filtering, upsampling, and signal delaying. 20 
These elements also perform interpolation of input state 
signals. For signal paths 2112 and 2114, this interpola- 
tion results in output signals 218 and 220 having the 
same clock rate. 

[0069] Position signal path 2114 includes a preproc- 25 
essing component 2140, a pre-interpolation filter 2142, 
an interpolation component 2144, and atime delay com- 
ponent 21 46. As described above, position signal path 
2114 receives trajectory data in the form of state signals 
210, 212, and 214, and generates position output signal so 
218. 

[0070] Pre-processing component 21 40 receives po- 
sition state signal 21 0, velocity state signal 212, and ac- 
celeration state signal 214. These signals are proc- 
essed to generate an integrated position signal 2150, 35 
which is at the slow dock rate. Integrated position signal 
2150 Is based on acceleration state signal 214, velocity 
state signal 212 and position state signal 210. To gen- 
erate integrated position signal 2150, preprocessing 
component 2140 performs a discrete time Integration 40 
operation, 

P„=P,+ nT*[l/)+(rjr*>A/2)l 

45 

As described above, P„ represents Integrated position 
signal 2150, P/ represents position state signal 210, V/) 
represents an velocity state signal 212, n is an integer 
representing a sample time at the slow clock rate, T is 
the slow clock period, and A represents acceleration so 
state signal 214. 

[0071] Pre-interpolation filter 2142 receives position 
state signal 210, and integrated position signal 2150. 
Pre-interpolation filter 21 42 processes these signals to 
generate filtered position signal 21 52. By generating fil- 55 
tered position signal 2152, pre-interpolation filter 21 42 
filters integrated position signal 2150 to limit the fre- 
quency content of signal 2152 before It is sent to inter- 



polation component 2144. Pre-interpolation filter 2140 
employs FIR filtering techniques that run at the slow 
clock rate. However, pre-interpolation filter 2140 can 
employ other types of discrete time filtering techniques, 
such as MR filtering. 

[0072] Interpolation component to 2144 receives fil- 
tered position signal 21 52 and position state signal 210. 
Upon receipt of these signals, Interpolation component 
21 44 perfomns two functions. First, interpolation compo- 
nent 2144 performs an upsampling operation to create 
a signal that is at the fast clock rate. This upsampling 
operation involves inserting L-1 zeros between each 
sample where L is the interpolation rate. 
[0073] Second, interpolation component 2144 per- 
forms a filtering operation on the upsampled signal. This 
filtering operation prevents high frequency resonance 
components from being Introduced by the upsampling 
function. The filtering Is perfomied using a low-pass FIR 
filtering operation. However, other types of filtering, such 
as IIRfilteringcan be employed. Together, these upsam- 
pling and filtering operations are collectively referred to 
herein as interpolation. The result of this Interpolation 
operation results in an interpolated position signal 21 54, 
which is sent to time delay component 2146. 
[0074] Both pre-interpolation filter 2142 and interpo- 
lation component 2144 employ "delta processing" tech- 
niques to increase the precision of the filtering and in- 
terpolation operations that they perfonn. These delta 
processing techniques are described below in greater 
detail with reference to FIGS. 22 and 23. 
[0075] Time delay component 2146 receives interpo- 
lated position signal 2154. This component operates at 
the fast clock rate and enables the propagation time and 
phase of position signal path 2114 to be adjusted. Time 
delay component 2146 generates position signal 218 
and sends it to output signal interface 2104. 
[0076] The description now turns to acceleration sig- 
nal path 2112. Acceleration signal path 2112 includes 
an interpolation component 21 20, and a time delay com- 
ponent 21 22. Acceleration signal path 21 1 2 receives ac- 
celeration state signal 214 from state signal Interface 
2102 and generates acceleration output signal 220. 
[0077] Interpolation component 2120, like interpola- 
tion component 2144, perfomris two functions. First, in- 
terpolation component 21 20 upsamples signal acceler- 
ation state signal 214 at an Interpolation rate of L. Sec- 
ond, interpolation component 2120 performs a filtering 
operation on this upsampled data. This filtering opera- 
tion prevents high frequency resonance components 
from being introduced by the upsampling function. As a 
result of these functions, interpolation component 2120 
produces an interpolated acceleration signal 2130 that 
is transferred to time delay component 2122. 
[0078] Time delay component 2122 receives interpo- 
lated acceleration signal 2130 and enables the propa- 
gation time and phase of acceleration signal path 2112 
to be adjusted. By incurring such delays, time delay 
component 2122 generates acceleration output signal 
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220, which is sent to output Interface 2104. 
Delta Processing 

[0079] Elements of profile executor 1 50, such as filter s 
2142 and Interpolation component 2144 employ a delta 
processing feature. This feature reduces computational 
errors in output signals, such as position output signal 
218. Delta processing involves subtracting position 
state signal 210 from time varying signals that represent io 
an absolute position. This subtraction results in "delta 
signals" that have a smaller dynamic range than the cor- 
responding signals that represent an absolute position. 
[0080] FIG. 22 is a block diagram illustrating an Im- 
plementation of filter 2142 that employs delta process- ^5 
Ing. This implementation includes an absolute position 
array 2202, an FIR coefficient array 2204, a first scalar 
summation node 2206, a vector product node 2208, a 
vector summation node 2210, a second scalar summa- 
tion node 2212, and a delay module 2214. 
[0081 ] Absolute position array 2202 receives an input 
signal (i.e., Integrated position signal 2150) from pre- 
processing component 2140. Absolute position array 
2202 processes the integrated position signal to provide 
a parallel position signal set 2230. Parallel position sig- 
nal set 2230 includes a plurality of discrete time signals. 
For each of the discrete time steps, these signals in- 
clude a plurality of contiguous values selected from sig- 
nal 2150. Thus, through parallel position signal set 
2230, absolute position array 2202 concurrently sends 30 
scalar summation node 2106 a plurality of contiguous 
values from signal 2150. 

[0082] Scalar summation node 2206 receives parallel 
position signal set 2230 and position state signal 21 0. 
Scalar summation node 2206 subtracts position state 35 
signal 21 0 from each of the signals within parallel signal 
set 2230. This subtraction produces a delta parallel sig- 
nal set 2232, which is sent to vector product node 2208. 
Signal set 2232 has a delta magnitude scale. This scale 
is smaller than the absolute magnitude scale of signal 
set 2230. 

[0083] Vector product node 2208 performs an ele- 
ment-by-element multiplication on delta parallel signal 
set 2232 and a plurality of FIR filter coefficients that are 
stored In FIR coefficient array 2204. This multiplication ^5 
produces a parallel product signal set 2234, which Is 
sent to vector summation node 2210. 
[0084] Vector summation node 221 0 sums each of the 
signals within parallel product signal set 2234. This sum- 
mation results in a single (Le. , a scalar) delta filtered sig- so 
nal 2236. 

[0085] Scalar summation node 221 2 receives delta fil- 
tered signal 2236 and a delayed position state signal 
2238. Delayed position state signal 2238 is position 
state signal 210 that is delayed a predetermined delay ss 
time by delay module 2214. This predetermined delay 
time matches the propagation delays associated with 
the above described transfer of signals between sum- 



mations node 2206 and 221 2. 

[0086] Scalar summation node 2212 adds delta fil- 
tered signal 2236 and delayed position state signal 
2238, thereby producing filtered position signal 2152. 
This addition returns filtered signal 2236 from a delta 
scale to an absolute scale. As described herein, filtered 
position signal 2152 is sent to interpolation component 
2144. 

[0087] The delta processing approach of this Imple- 
mentation enables filter 2142 to minimize round up er- 
rors incurred by computations (I.e., summations and 
multiplications) that are associated with filtering opera- 
tions. This en^or minimization is achieved through the 
subtraction of position state signal 21 0 from parallel po- 
sition signal set 2230 at scalar summation node 2206. 
This subtraction results in signal set 2232, which has a 
smaller magnitude than signal set 2230. 
[0088] In fixed precision processing environments, 
round up error due to fixed precision numeric represen- 
tations are a direct function of the magnitude of input 
signals. By reducing the magnitude of input signals, as 
in delta parallel signal set 2232, the computations per- 
formed at nodes 2208 and 2210 are minimized. This 
minimization yields an output signal (i.e., filtered posi- 
tion signal 2152) that has smaller errors. 
[0089] FIG. 23 is a block diagram illustrating an im- 
plementation 500 of interpolation component 21 44. This 
Implementation is similar to the above Implementation. 
However, this implementation includes an upsampling 
module 2302 that increases the sampling rate of an in- 
put signal (I.e., filtered position signal 2152), thereby 
producing an upsampled signal 2322, which is sent to 
absolute position anray 2303. This implementation also 
receives position state signal 21 0 in order to employ the 
delta processing feature that is described above with 
reference to FIG. 22. Thus, as in FIG: 21, the delta 
processing technique employed by this implementation 
component produces an output signal (i.e., interpolated 
signal 2154) that has smaller errors. 
[0090] Thus, the present invention includes a tech- 
nique for efficiently processing a discrete time input sig- 
nal, p(n), having a first clock rate Into a discrete time 
output signal having a second clock rate. This technique 
minimizes computational errors. 

111. Conclusion 

[0091] While various embodiments of the present In- 
vention have been described above, It should be under- 
stood that they have been presented by way of example, 
and not limitation. It will be apparent to persons skilled 
in the relevant art that various changes In detail can be 
made therein without departing from the spirit and scope 
of the invention. Thus the present invention should not 
be limited by any of the above-described exemplary em- 
bodiments. 
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Claims 

1 . A method of controlling the movement of a compo- 
nent of a servo-driven photolithography tool, com- 
prising: 5 

(a) receiving initial control data regarding re- 
quired start and end points, distances, and ve- 
locities of scans; 

(b) planning an infinite jerk, multi-axis trajectory io 
of the component on the basis of the initial con- 
trol data; 

(c) perfonning, in real time, a profile execution 
on the basis of the planned trajectory, wherein 

the trajectory is shaped using digital filtering; i5 
and 

(d) outputting, to a control system, execution 
data that comprises positions and accelera- 
tions of the component. 

20 

2. The method of claim 1 , wherein the component is a 
wafer stage. 

3. The method of claim 1 , wherein the component is a 
reticle stage. 2S 

4. The method of claim 1 , wherein said step (b) com- 
prises: 

(i) calculating a preliminary motion trajectory for 30 
each axis; 

(li) adjusting zero velocity time intervals in the 
preliminary motion trajectories, if necessary, to 
synchronize motion across axes; and 
(ill) inserting delays at the start of the prelimi- 35 
nary motion trajectory for each axis to compen- 
sate for fitter delays. 

5. The method of claim 4, wherein said step (i) com- 
prises: 40 

(A) adjusting a constant velocity scan length to 
account for filter width, settling time, and quan- 
tization; 

(B) computing acceleration intervals between 
consecutive constant velocity scans; 

(C) creating a zero velocity point between the 
consecutive constant velocity scans; and 

(D) quantizing the acceleration intervals. 

6. The method of claim 5, wherein said step (B) com- 
prises: 

calculating an optimal velocity profile for each 

axis; 

If maximum velocity or maximum acceleration 
limits are exceeded, constraining the optimal 
velocity profile according to the limits; and 



calculating transition times between consecu- 
tive constant velocity scans. 

The method of claim 5. wherein said step (C) com- 
prises: 

if a trajectory between scans either starts or 
ends at zero velocity, or has no reversal point, 
brealcing the trajectory into two sub-optimal 
segments, where the first sub-optimal segment 
ends with zero velocity and the second sub-op- 
timal segment begins with zero velocity; 
otherwise, calculating a reversal point and split- 
ting the trajectory at the reversal point, into two 
segments, wherein velocity is zero at the re- 
versal point. 

8. The method of claim 4, wherein said step (ii) com- 
prises: 

(A) identifying a non-critical axis trajectory; 

(B) calculating a time difference between a crit- 
ical axis trajectory and the non-critical axis tra- 
jectory; and 

(C) introducing a delay, equal to the time differ- 
ence, at a zero velocity point of the non-critical 
axis trajectory. 

9. The method of claim 4, wherein said step (b) further 
comprises the following step, perfomied after said 
step (ii): 

(iv) inserting a delay at the start of an axis tra- 
jectory to compensate for filter delays. 

10. The method of claim 9, wherein said step (iv) com- 
prises: 

(A) identifying an axis trajectory with the long- 
est execution filter delay; 

(B) identifying an axis trajectory with a shorter 
execution filter delay; 

(C) calculating the difference between the exe- 
cution filter delays; and 

(D) introducing a matching delay, equal to the 
difference between the execution filter delays, 
into the axis trajectory with the shorter execu- 
tion fitter delay. 



50 11. The method of claim 1 , wherein said step (C) com- 
prises: 

(i) receiving trajectory data; 

(ii) creating an integrated position signal based 
55 on the trajectory data; 

(iii) filtering the integrated position signal, using 
an initial position signal from the trajectory data 
to create a filtered position signal; 
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(iv) interpolating the filtered position signal, us- 
ing the initial position signal, to create an inter- 
polated position signal; and 
(V) delaying the interpolated position signal to 
compensate for propagation delay and phase s 
shift; 

(vl) interpolating an acceleration state signal 
from the trajectory data to create an Interpolat- 
ed acceleration signal; and 
(vii) delaying the interpolated acceleration sig- io 
nal to compensate for propagation delay and 
phase shift. 

12. The method of claim 1 1 , wherein said step (iii) com- 
prises: 

(A) applying a first absolute position array to the 
initial position signal to create a parallel signal 
set; 

(B) subtracting the initial position signal from 
each parallel position signal of the parallel po- 
sition signal set to create a delta parallel posi- 
tion signal set; 

(C) multiplying the delta parallel position signal 
set by first filter coefficients to create a parallel 
product signal set; 

(D) summing all parallel product signals of the 
parallel product signal set, to create a scaler 
delta filtered signal; and 

(E) adding a first delayed initial position signal 30 
to the delta filtered signal to create a filtered po- 
sition signal. 

13. The method of claim 1 1 , wherein said step (iv) com- 
prises: 35 



the interpolated signal. 

15. The method of claim 1 1 . wherein said step (vl) com- 
prises: 

(A) upsampling the acceleration state signal; 
and 

(B) filtering the upsampled acceleration signal. 

16. A system for controlling movement of a component 
of a servo driven photolithography tool, comprising: 

a trajectory planning module that plans an infi- 
nite jerk, multi-axis trajectory for the compo- 
nent; and 

a profile executor module that plans execution 
of said trajectory in a manner that optimizes 
said trajectory, wherein said executor module 
comprises a digital filter. 

1 7. The system of claim 1 6, wherein said executor mod- 
ule comprises: 

a first interpolation component that receives an 
acceleration signal from said trajectory plan- 
ning module and Interpolates said acceleration 
signal; 

a preprocessor that receives said acceleration 
signal, and velocity and initial position signals 
from said trajectory planning module and cre- 
ates an integrated position signal; 
a filter that smoothes said integrated position 
signal to create a filtered position signal; and 
a second interpolation component that receives 
and interpolates said filtered position signal. 



(A) upsampling the filtered position signal; and 

(B) filtering the upsampled filtered position sig- 
nal to create the interpolated signal. 

40 

14. The method of claim 1 3, wherein said step (B) com- 
prises: 

applying a second absolute position array to the 
upsampled filtered position signal to create a 
parallel upsampled position signal set; 
subtracting the initial position signal from each 
parallel upsampled position signal of the paral- 
lel upsampled position signal set to create a 
delta parallel upsampled position signal set; so 
multiplying the delta parallel upsampled posi- 
tion signal set by second filter coefficients to 
create a parallel product upsampled signal set; 
summing all signals of the parallel product up- 
sampled signal set to create a scalar delta fll- 55 
tared up sampled signal; and 
adding a second delayed Initial position signal 
to the delta filtered upsampled signal to create 



18. A computer program product comprising a compu- 
ter useable medium having control logic stored 
therein for causing a computer to control the move- 
ment of a component of a servo-driven photolithog- 
raphy tool, the computer control logic comprising: 

(a) first computer readable program code 
means for causing the computer to receive ini- 
tial control data regarding required start and 
end points, distances, and velocities of scans; 

(b) second computer readable program code 
means for causing the computer to plan an in- 
finite jerk, multi-axis trajectory of the compo- 
nent on the basis of the initial control data; 

(c) third computer readable program code 
means for causing the computer to perform, in 
real time, a profile execution on the basis of the 
planned trajectory, wherein the trajectory is 
shaped using digital filtering; and 

(d) fourth computer readable program code 
means for causing the computer to output exe- 
cution data that comprises positions and accel- 
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erations of the component, to a control system. 

1 9. The computer program product of claim 1 8, wherein 
the component Is a wafer stage. 

5 

20. The computer program product of claim 1 8. wherein 
the component is a reticle stage. 

21 . The computer program product of claim 1 8, wherein 
said second computer readable program code io 
means comprises: 

(i) computer readable program code means for 
causing the computer to calculate a preliminary 
motion trajectory for each axis; is 
(il) computer readable program code means for 
causing the computer to adjust zero velocity 
time Intervals, In the preliminary motion trajec- 
tories, If necessary, to synchronize motion 
across axes; and 20 
(lii) computer readable program code means 
for causing the computer to insert delays at the 
start of the preliminary motion trajectory for 
each axis, to compensate for filter delays. 

25 

22. The computer program product of claim 21 , wherein 
said computer readable program code means (i) 
comprises: 

(A) computer readable program code means 30 
for causing the computer to adjust constant ve- 
locity scan length to account for filter width, set- 
tling time, and quantization; 

(B) computer readable program code means 

for causing the computer to compute accelera- 35 
tion intervals between consecutive constant ve- 
locity scans; 

(C) computer readable program code means 
for causing the computer create a zero velocity 
point between the consecutive constant veloc- ^0 
Ity scans; and 

(D) computer readable program code means 
for causing the computer to quantize the accel- 
eration intervals. 

45 

23. The computer program product of claim 22, wherein 
said computer readable program code means (B) 
comprises: 

computer readable program code means for so 
causing the computer to calculate an optimal 
velocity profile for each axis; 
computer readable program code means for 
causing the computer to constrain the optimal 
velocity profile according to the limits if maxi- 55 
mum velocity or maximum acceleration limits 
are exceeded; and 

computer readable program code means for 



causing the computer to calculate transition 
times between consecutive constant velocity 
scans. 

24. The computer program product of claim 22, wherein 
said computer readable program code means (C) 
comprises: 

computer readable program code means for 
causing the computerto break the trajectory in- 
to two sub-optimal segments, if a trajectory be- 
tween scans either starts or ends at zero veloc- 
ity, or has no reversal point, where the first sub- 
optimal segment ends with zero velocity and 
the second sub-optimal segment begins with 
zero velocity; 

computer readable program code means for 
causing the computer to calculate a reversal 
point and splitting the trajectory at the reversal 
point, into two segments otherwise , wherein ve- 
locity is zero at the reversal point. 

25. The computer program product of claim 21 , wherein 
said computer readable program code means (ii) 
comprises: 

(A) computer readable program code means 
for causing the computer to identify a non-crit- 
ical axis trajectory; 

(B) computer readable program code means 
for causing the computer to calculate a time dif- 
ference between a critical axis trajectory and 
the non-critical axis trajectory; and 

(C) computer readable program code means 
for causing the computer to introduce a delay, 
equal to the time difference, at a zero velocity 
point of the non-critical axis trajectory. 

28. The computer program product of claim 21 , wherein 
said second computer readable program code 
means further comprises: 

(iv) computer readable program code means 
for causing the computerto insert a delay at the 
start of an axis trajectory, to compensate for fil- 
ter delays. 

27. The computer program product of claim 26, wherein 
said computer readable program code means (iv) 
comprises: 

(A) computer readable program code means 
for causing the computer to identify an axis tra- 
jectory with the longest execution filter delay; 

(B) computer readable program code means 
for causing the computer to identify an axis tra- 
jectory with a shorter execution filter delay; 

(C) computer readable program code means 
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for causing the computer to calculate the differ- 
ence between the execution filter delays; and 
(0) computer readable program code means 
for causing the computer to introduce a match- 
ing delay, equal to the difference between the s 
execution filter delays, into the axis trajectory 
with the shorter execution filter delay. 



for causing the computer to sum all parallel 
product signals of the parallel product signal 
set, to create a scaler delta filtered signal; and 
(E) computer readable program code means 
for causing the computer to add a first delayed 
initial position signal to the delta filtered signal, 
to create a filtered position signal. 



28. The computer program product of claim 1 7, wherein 
said computer readable program code means (C) 

comprises: 

(i) computer readable program code means for 
causingthecomputerto receive trajectory data; 

(ii) computer readable program code means for 
causing the computer to create an integrated 
position signal based on the trajectory data; 
(III) computer readable program code means 
for causing the computer to filter the intergrated 
position signal, using an initial position signal 
from the trajectory data, to create a filtered po- 
sition signal; 

(iv) computer readable program code means 
for causing the computer to interpolate the fil- 
tered position signal, using the initial position 
signal, to create an interpolated position signal; 

(v) computer readable program code means for 
causing the computer to delay the interpolated 
position signal to compensate for causing the 
computer propagation delay and phase shift; 

(vi) computer readable program code means 
for causing the computer to Interpolate an ac- 
celeration state signal from the trajectory data, 
to create an interpolated acceleration signal; 
and 

(vii) computer readable program code means 
for causing the computer to delay the interpo- 
lated acceleration signal, to compensate for 
propagation delay and phase shift. 

29. The computer program product of claim 28, wherein 
said computer readable program code means (iii*) 
comprises: 

(A) computer readable program code means 
for causing the computer to apply a first abso- 
lute position array to the initial position signal, 
to create a parallel signal set; 

(B) computer readable program code means 
for causing the computer to subtract the initial 
position signal from each parallel position sig- 
nal of the parallel position signal set. to create 
a delta parallel position signal set; 

(C) computer readable program code means 
for causing the computer to multiply the delta 
parallel position signal set by first filter coeffi- 
cients, to create a parallel product signal set; 

(D) computer readable program code means 



30. The computer program product of claim 28, wherein 
10 said computer readable program code means (Iv) 

comprises: 

(A) computer readable program code means 
for causing the computer to upsample the iW- 

IS tered position signal; and 

(B) computer readable program code means 
for causing the computer to filter the up sam- 
pled filtered position signal, to create the inter- 
polated signal. 

20 

31 . The computer program product of claim 30, wherein 
said computer readable program code means (B) 
comprises: 

25 computer readable program code means for 

causing the computer to apply a second abso- 
lute position an^ay to the up sampled filtered po- 
sition signal, to create a parallel up sampled po- 
sition signal set; 
30 computer readable program code means for 

causing the computer to subtract the Initial po- 
sition signal from each parallel up sampled po- 
sition signal of the parallel up sampled position 
signal set, to create a delta parallel up sampled 
35 position signal set; 

computer readable program code means for 
causing the computer to multiply the delta par- 
allel up sampled position signal set by second 
filter coefficients, to create a parallel product up 
40 sampled signal set; 

computer readable program code means for 
causing the computer to sum all signals of the 
parallel product up sampled signal set, to cre- 
ate a scalar delta filtered up sampled signal; 
43 and 

computer readable program code means for 
causing the computer to add a second delayed 
Initial position signal to the delta filtered up 
sampled signal, to create the interpolated sig- 
50 rial. 

32. The computer program product of claim 28, wherein 
said computer readable program code means (vi) 
comprises: 

55 

(A) computer readable program code means 
for causing the computer to upsample the ac- 
celeration state signal; and 
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(B) computer readable program code means 
forcausing the computer to filter the upsampled 
acceleration signal. 
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